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The electromagnetic interaction with quarks is investigated through a relativis- 
tic, electromagnetic gauge-invariant treatment. Gluon dressing of the quark-photon 
vertex and the quark self-energy functions is described by the inhomogeneous Bethe- 
Salpeter equation in the ladder approximation and the Schwinger-Dyson equation 
in the rainbow approximation respectively. Results for the calculation of the quark- 
photon vertex are presented in both the time-like and space-like regions of photon 
momentum squared, however emphasis is placed on the space-like region relevant 
to electron scattering. The treatment presented here simultaneously addresses the 
role of dynamically generated qq vector bound states and the approach to asymp- 
totic behavior. The resulting description is therefore applicable over the entire range 
of momentum transfers available in electron scattering experiments. Input param- 
eters are limited to the model gluon two-point function, which is chosen to reflect 
confinement and asymptotic freedom, and are largely constrained by the obtained 
bound-state spectrum. 



I. INTRODUCTION 

As hadrons are subjected to more detailed examination in electron scattering experi- 
ments, such as those proposed at CEBAF, a quark-based description of their electromagnetic 
(EM) interaction, which is applicable over a broad range of momenta, becomes desirable. 
Local effective field theories based on hadronic degrees of freedom, although attractive for 
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their efficiency, are subject to difficulties when conditions are sufficient to probe charac- 
teristics of subhadronic origin M. An ideal perspective on this problem can be gained by 
sacrificing local interactions in favor of retaining some knowledge of the quark substructure 
content of the effective hadron fields 0, and is the spirit of the present investigation. In such 
a treatment, the EM interaction with a pion at tree level in the hadron fields is, for example, 
illustrated in Fig.|l[ There, the pion Bethe-Salpeter amplitudes, T, reflect the ability of the 
internal degrees of freedom to share the transferred momentum, while the quark-photon 
vertex (QPV), T u , contains information about the interaction with a dressed quark, includ- 
ing the production of qq vector-meson modes expected from vector-meson dominance or 
dispersion theory phenomenology. At larger momentum transfers where the effective scale 
provided by the vector-meson poles has been exceeded, the QPV also provides information 
about the approach of photon-hadron vertices to asymptotic behavior. A quantitative un- 
derstanding of the QPV is therefore central to the EM description of hadrons at intermediate 
momentum scales, and must address issues of confinement, dynamically generated qq vector 
bound states, and asymptotic freedom. 

In addition to its role in the EM interaction, the QPV is a quantity of interest on its 
own merit. It is well known, for example, that a formal description of meson bound-state 
spectra and wavefunctions can be obtained by considering the interaction of a quark with 
an external field having the quantum numbers of the bound state of interest. In such a 
description, the external field excites a qq pair from the vacuum. Their subsequent gluon- 
mediated interactions provide self-energy and vertex dressing. The existence of bound states 
is exhibited by poles in the dressed vertex. The position and residue of a given pole yield 
information about the bound-state mass and wavefunction respectively. Vertex functions 
of quarks with external fields thus provide a unique environment in which to study the 
formation of meson bound states. The interaction with a photon is simply one example of 
this, and it is expected that the approach presented here can be applied to the description 
of bound states with other quantum numbers as well as vectors. 

In this paper the QPV with gluon dressing at the ladder level of approximation is inves- 
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tigated in both the time-like and space-like region of photon momentum squared. Particular 
attention is given to the formation of vector bound states and confinement phenomenol- 
ogy, however the primary goal of this initial investigation is to develop a description of the 
quark-photon interaction in the space-like region relevant to electron scattering which si- 
multaneously addresses the role of vector mesons and quark substructure. The equations 
implied by the ladder approximation are 1 the inhomogeneous Bethe-Salpeter equation 

r„(P, q) = -i ltl - ^g 2 J 0^D{P - K) lv G{K + )Y„{K, q)G{KJ) lu , (1) 

and the Schwinger-Dyson equation 

s(p) = t g2 1 Wf D{p " K ^ K ^ w 

where K± = K ± g/2, the quark Green's function, G, is defined from its inverse 
G^ 1 (P) = iP + Y,(P), and for convenience the model gluon two-point function is taken 
to be diagonal in Lorentz indices (D^ = S^D). In Eq.(|I]) q is the photon momentum, 
while P is the average of the incoming and outgoing quark momenta. These expressions are 
illustrated in Fig.|[ The current-quark masses in this investigation are assumed to be zero, 
which allows an overall charge matrix to be omitted from Eq. ([!]), but the approach can be 
easily extended to include nonzero masses. The absence of masses further implies that all 
dimensionful constants enter through the parameterization of the gluon two-point function, 
D. Here the infrared contribution to this quantity is characterized by an effective strength, 
rj, and an effective range, Rq, in addition to the known ultraviolet form. 

The motivation for employing the ladder approximation is provided by the following 
arguments. First, it offers an EM gauge- invariant description. Both the Ward-Takahashi 
identity (WTI), 

q fl T fl (P,q) = G- 1 (P_)-G- 1 (P + ), (3) 
and the Ward identity (WI), 
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are satisfied by this approximation. This is easily demonstrated by directly substituting 
(0) or (JD into (|1]) for the vertex, and employing (0). Further, it has been shown recently 
how gauge invariance at the quark level in the ladder approximation to the QPV fits into a 
gauge-invariant description of composite mesons leading to WTIs at the hadronic level 0. 
Second, in the vicinity of a qq resonance, the vertex obtains the form 

M 2 

where is the residue at the pole, and from (|l|) satisfies 

ty,(P, q) = -\g 2 J 0^ D ( P - K) lv G{K+)n»{K, q)G{KJ) lv , (6) 

for q 2 = —My. This is the homogeneous Bethe-Salpeter equation in the ladder approxima- 
tion for a vector qq bound state. One can conclude from (j3|) and (|5D that the on-mass-shell 
solution of the homogeneous equation satisfies g^f^P, g) = 0, provided the self-energy 
functions are nonsingular. The Bethe-Salpeter equation at this level has been investigated 
previously, and acceptable results have been obtained for the lowest mass bound states ||-|| . 
It is therefore anticipated that, in the addition to providing an EM gauge invariant descrip- 
tion of the quark-photon interaction for use in electron scattering, this level of approximation 
has the ability to describe the spectrum and structure of vector bound states. 

Finally, the Schwinger-Dyson equation has been extensively studied with regard 
to confinement phenomenology and dynamical chiral symmetry breaking 0. It has been 
demonstrated, for example, that for a suitably chosen gluon two-point function D, quark 
Green's functions which are free of singularities on the real P 2 axis are produced; therefore 
denying the existence of a single-particle spectrum. It is shown here, however, that this is not 
a sufficient condition for the formation of finite-size qq states, but is an essential ingredient in 
a description of confinement based on an interaction between dressed quarks which does not 
rise to infinity at large distances. This alternate description of confinement in fact relies on 
the interaction between dressed quarks diminishing at separations beyond the characteristic 
length scale, Ro, at which point the quark Green's function approaches its vacuum value. 
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The quarks are then repelled by the vacuum due to absence of a mass pole in the vacuum 
quark Green's function. The effective "rising potential" of traditional models of confinement 
occurs here as a result of the repulsive interaction of quarks with the vacuum in combination 
with the decreasing interaction between dressed quarks with increasing separation. Such a 
description has proven successful in the formulation of a confining nontopological soliton 
model for baryons |§. This mechanism arises naturally here, and is pursued as a possible 
realization of confinement. 

The paper is organized as follows. In Sec. 2 the general features of the QPV are given 
along with the reduction of Eq.([l]) to a form suitable for computer applications. The so- 
lution of the homogeneous Bethe-Salpeter equation (H) is also discussed there. In Sec. 3 a 
discussion of confinement is presented in the context of a particular model which yields an- 
alytic solutions to both the inhomogeneous Bethe-Salpeter and Schwinger- Dyson equations 
(H) and (Q) respectively. Also discussed there is the criteria for the formation of finite-range 
bound states in the context of this model of confinement. The numerical evaluation of the 
QPV for both the time-like and space-like regions of the photon momentum squared is given 
in Sec. 4. Finally, a summary is presented in Sec. 5. 



In this section, the steps required to reduce the inhomogeneous Bethe-Salpeter equation 
(|I|) to a form suitable for computer applications are outlined. The approach is straightfor- 



lengthy algebraic manipulations are omitted. To begin, it is useful to write the self energy, 



II. GENERAL FEATURES 



A. The inhomogeneous Bethe-Salpeter equation 



ward and is based on the general features of the QPV as described by ([T]). Where possible, 





(7) 



where 
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>2 -8 f d*K n/D _ A(K 2 )P-K 



A(P 2 ) - 1 P 2 = g 2 - - — -D(P - AT) - o ; ox ; — — — , (8) 

V ' J y 3 7 (2vr) 4 V ! K 2 A 2 {K 2 ) + B 2 {K 2 Y KJ 

and 

PfP 2 ) - a 216 ( * K D(P K) B{R2) (9) 

B{P) - 9 YJ { ~ ] K 2 A 2 (K 2 ) + B 2 (K 2 Y (9) 

The quark Green's function, G, and its inverse, can then be written as 

G(P) = -i Pa{P 2 ) + f3{P 2 ) and G~\P) = i PA(P 2 ) + B(P 2 ), (10) 

where a and (3 are defined in an obvious fashion from G^ 1 . 

From Eq. (|IJ) it is evident that the QP V, in addition to its Lorentz four- vector structure, 
is a four-by-four matrix in Dirac indices and is in general a function of two non-orthogonal 
four vectors, P and q (P ■ q ^ 0). On first appraisal, one might therefore expect that 
the complete set of sixteen Dirac matrices {1,75,7^,757^,0^} must be employed in its 
description. However, by limiting the gluon two-point function, D, to diagonal components 
in Lorentz indices, the tensor, a^, is excluded. This is seen explicitly by rewriting the 
integrand of ([!]) using 

[ lv G{K + )Y,{P,q)G{KJ) lv ] u = ( lw ) ij [G(K + )T,(P,q)G(K_)] jk ( lu ) kl 

= ij^gtr [M a G(K + )T,(P, q)G(K_)} , (11) 

where the set of matrices M a = jl, 175, 47^/^/2, i757 At /\/2|, from which a^ u is absent, are 
obtained through Fierz reordering, and the repeated indices are summed. Further, the 
pseudo-scalar matrix, 75, can be excluded by considering that the corresponding term in 
(|TT|) has the form 7sA^(P, q), where the matrix structure has been completely factored. In 
order to contribute to the QPV, this term must form a four vector, which requires A^(P, q) 
to be an axial vector. It is immediately evident that this is impossible since there is an 
insufficient number of four vectors (P and q) to combine with the tensor to form an 

axial vector. The solution to ([I]) must therefore have the general form 

T,(P, q) = 1A«(P, q) + z T ,Ag(P, q) + i lblu e^P a q^ 2 A^ (p, q ), (12) 
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where the matrix structure is now explicit and the quantities Aw are dimensionless functions 
with Lorentz structure indicated by their indices. The constant rj has dimension of (length) _1 
and is associated with the infrared strength of the gluon two-point function. 

The further reduction of the AW to a set of invariant functions is achieved through the 
use of the symmetry transformations of the vertex r M (P, q) under 75; charge conjugation, 
C = 7274; and parity, V = 74. The transformation properties are determined directly from 
the inhomogeneous Bethe-Salpeter equation (|1|) and are given by 

7 5 r M (P,g)7 5 = -r At (-P,-g) 

cryp^c-^-ry-pg)' (13) 
pryp, q )v~ l = ^r i/ (p 4 ; -p, g 4 ; -q) 

where uo^ u =diag(— 1, —1, —1, 1) and t denotes a matrix transpose. From the general form 
given in (|l^) and the first two expressions in flHf), it is clear that Affi(P,q) and A®(P, q) 
are even in both P and g, while A^(P, q) is odd in P and even in q. The quantities AW can 
therefore be written as 

AW(P,g) = &^Af + 3: A T 

AS(P,*) = ^^ + ^AJ (14) 



A (3) (P,g) = A 



Here the coefficients, A, are functions of P 2 , g 2 , and Cp where Cp q = P • q/(Pq) is the 
direction cosine between P and g, and Pj = P M — P • gg^/g 2 is the vector transverse to 
g^ (P T • g = 0). The behavior of the coefficients for large space-like g 2 consistent with 
asymptotic freedom |9j can be anticipated from ( |T2D and (0). In particular, as g 2 —>■ oo the 
coefficients X^' L must tend to unity, while all others tend to zero. The separation in flU]) 
into longitudinal and transverse contributions with respect to the index ^ allows three of 
the eight coefficient functions, A, to be immediately determined from the WTI @ in terms 



of the self-energy functions A and B given in 
coefficient functions are thus given by 



and (pi) respectively. The longitudinal 



P ' q -\[ = B(P 2 ) - B(Pl) 



V 



(15) 



A 



A{P 2 _) + A{P_ 



The additional constraint provided by the WI (||) obtains the soft-photon limits: 





9=0 ~ 


Af 


g=0 


= -2r)B'(P 2 ) 


X L 2 


9=0 " 


A2 


9=0 


= -2rfA'{P 2 


Af 


9=0 ~ 


Af 


9=0 


= A(P 2 ), 



(16) 



where the prime denotes differentiation with respect to the argument. 



The five unknown transverse coefficients can be isolated using ([12]) and (14), and the 
orthogonality of the Dirac matrices. The closed set of equations so obtained can be written 
in matrix form as 

b^MAP^c^-g^xJiP 2 ,^^ 2 ) 

+ j Q dKK 3 dC Kq M l3 (P 2 , C 2 Pq , K 2 , C 2 Kq , q 2 )\j(K 2 , C 2 Kq , q 2 ), (17) 

where Ai and Af(b and A T ) are matrices (vectors) in the five- component space of the coef- 
ficient functions and Ck q is the direction cosine between K and q. The steps required to 
reach the result ([17]) along with the explicit forms of A4, M and b are given in the appendix. 
The equation (|I7]) contains the same information as the transverse component of the inho- 
mogeneous Bethe-Salpeter equation ([]]); to this point, no approximations have been made. 
To proceed further requires the specification of the gluon two-point function, D. 



B. The homogeneous Bethe-Salpeter equation 



The solutions of the homogeneous Bethe-Salpeter equation, (|6|), describe the coupling 
of a qq vector bound state to a quark, and thereby provide information about the internal 
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dynamics of the bound state. Although these solutions can be obtained directly from the 
coefficient functions, Xj of the previous section, their normalization is more readily achieved 
by considering the solution of the homogeneous equation. To this end, the quantity f2 M (P, q) 
is decomposed as 

n„(p, q ) = i$\p, q) + Af^(P, q) + A^ibi^(p, q ). (is) 

The association of the coefficients (fffl with the coefficients of (0) is directly observed 
from the relationship between the solutions to the inhomogeneous and homogeneous Bethe- 
Salpeter equations given in (|5|). The substitution of flI8|) into the homogeneous equation (||) 
allows a description in the form of the eigenvalue equation 

J d 4 KA- b \P,K;q)^ b (K,q) = a(g 2 )<^(P,g) (19) 

to be obtained. The quantity A~ b (P, K; q) is defined as 

A- ab \P, K;q) = \j ^y^^ &ab + 5 (4 \P - K)ti {M a G(K + )M b G(K_)\ (20) 

and plays the role of an inverse propagator for the qq composite system, as has been demon- 
strated in previous work on the hadronization of quark field-theory models [jn|. That 
description is extended here to include the more general structure in QTSfl. The eigenfunc- 
tions are defined in terms of the coefficient functions, 0W, as $^ = <p^\ $2-5 = ^avi anc ^ 
$g_ g = . The matrices M a = l^ljr |l, 27^/ ilhln/ v^2~j include unity in the color and 
two-component flavor space, and the trace is taken to include their indices as well as those of 
the Dirac matrices. The homogeneous equation (|6|) is regained through the requirement that 
the eigenvalue vanishes on the mass shell, that is, a(q 2 = —M v ) = 0. The normalization 
condition for the eigenfunctions is obtained by multiplying ( |l"8D on the left by $%(P, —q) 
and integrating over relative momentum P, which gives 

J d'P^KZZiP, -q)A-}(P, K; q)& b (K, q) = a(q 2 ) J d l P<t>'(P, -q)<K(P, q) 

^(t„-!f\{t + l4)z{f). (21) 



The last line in ( pO|) follows from the fact that the on- mass-shell eigenfunctions, con- 
tain only components transverse to the center-of-mass momentum q^. The (dimensionless) 
wavefunction renormalization constant Z(q 2 ) can be absorbed into the Bethe-Salpeter am- 



plitudes by defining $^ = J Z(—M v ). The normalization condition can then be written 
explicitly as 



~ / *P#K& a {P, -q) A [A- b \P, K- q)] %{K, q) 



= 2q„ (22) 

-Ml 



which is the standard result JTT]]. In terms of the properly normalized Bethe-Salpeter am- 
plitude, the relationship to the QPV given in ([5]) now reads 

T(Pq)~ ^ (23) 



where fv = l/yZ(— M v ) is the effective coupling of the vector meson to the photon. 

III. A SIMPLE MODEL 

Traditional descriptions of quark confinement are based on a potential between quarks, 
in an appropriate color combination, which rises with increasing separation (to infinity in 
the absence of pair creation). Explored here is the notion that this effect may be in part 
due to the interaction of quarks (or more generally colored objects) with the vacuum. In 
particular, that confinement is manifest in the absence of a mass pole in the vacuum Green's 
function of a colored object implies that the object is repelled by the vacuum, and hence 
attracted to other colored objects to form a color singlet by virtue of the interaction with 
the vacuum. An illustration of this mechanism is afforded through the use of a simple model 
12, 8| of the gluon two-point function given by 

2^4 



D 



P _K)2) = ^!L!L 5 (4) ( p_ K) . (24) 



9 2 



At the present level of approximation, the self-energy dressing from the vacuum is provided 
by the equations (||) and (H). The solutions to these implied by (|4j) are 
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B(P 2 ) 



(jf -4P 2 )I, 




(25) 




That the solutions ( p5f) produce a model of confinement as described above can be seen 

by the absence of a solution to the expression P 2 + M 2 (P 2 ) = 0, with M = B/A. The 
presence of other colored objects in a color singlet configuration can supply the additional 
interaction necessary for the formation of a propagating mode, as is demonstrated in the 
nontopological soliton model of Ref. ||. There, for example, in the case of quarks coupled to 
a constant scalar mean field by the self-energy function B, the quark inverse Green's function 
obtains the form G^ 1 (P) = i PA(P 2 ) + B(P 2 )(1 + x), where > x > — 2 characterizes 
the strength of the mean field. In this case a continuous single-particle energy spectrum, 
E 2 (P 2 ) = P 2 + M 2 ., is obtained because of the infinite-range potential, with the constituent 
mass given by 



From (^) it is evident that the increase of the constituent mass as the strength, x, of the 

mean field decreases toward its vacuum value (x = 0), is due to the repulsive interaction with 
the vacuum. For the case of the constant mean field, the quarks are allowed to propagate 
throughout space. However, due to the energy stored in the mean field, the self-consistent 
(minimum-energy) solution acquires a finite range. In that case, as a quark in the system is 
separated from the others, the influence of the mean field on the quark diminishes and the 
mass rises as in (^6|). The quarks are thus confined to the region of nonzero mean field by 
virtue of their interaction with the vacuum. 

The situation is similar for the QPV of interest here. The use of the gluon two-point 
function ([24]) reduces the equations in ( |T7|) to algebraic form, which can then be solved by 
numerical or symbolic techniques. Shown in Fig. 3 is the solutions for the five transverse 




(26) 
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coefficient functions in both the time-like and space-like regions of photon momentum q 2 for 
P 2 = Cp q = 0. The single parameter, and the only dimensionful constant, rj, is given the 
value 5 fm _1 . It appears from Fig.3a that there is a pole for q 2 = —rf/2. Closer analysis of 
the q 2 — P 2 plane shows, however, that the singularity structure is not a discrete pole, but 
is instead a continuous spectrum. For the region P 2 < rf /4, the analytic solutions have the 
simple form 

Ar] \B(P 2 ) + B(P 2 ] 



7] 2 - AP 2 + 5q 2 + B(P 2 )B(P 2 

„2 



A 2 ^A 4 



AT — — 



2r] 2 - AP 2 + q 2 + B{Pl)B{P 2 ) 



A 



3r] 2 _ gp2 + Q g 2 

— 32?7 4 



(27) 



[r] 2 - AP 2 + 5q 2 + B(P 2 )B(Pl)} [3r] 2 - 8P 2 + 6q 2 } 

2 



5 3r) 2 - 8P 2 + Qq 2 

from which the singularity structure is visible. The self-energy function B is given in 
The analytic solutions outside this region are not as simple, but can be represented graph- 
ically. As an example, the coefficient function Aj is shown in Fig.|] for Cp q = 0. The 
exhibited continuous spectrum arises from the fact that the gluon two-point function (p4|) 
provides a constant, or infinite-range, interaction between the quarks in coordinate space. 
The quarks are therefore allowed to propagate throughout space uninhibited by the presence 
of the vacuum. This situation is analogous to that of the constant mean-field case discussed 
previously. 

Several observations can be made from this example. First, the absence of a mass pole in 
the vacuum Green's function of a colored object is sufficient for a description of confinement 
provided the subsequent interactions with other colored objects contain a scale beyond which 
the interaction diminishes allowing the presence of the vacuum to become apparent. This 
defines a confinement scale that is largely responsible for the characteristics of bound states. 
The effect of such a scale is investigated in the numerical work to follow in Sec. 4, however, 
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one can anticipate that as the range of the interaction increases, the density of bound-state 
poles increases; finally giving rise to a continuous spectrum as illustrated above. Second, 
the identification of a true bound-state pole requires careful examination of the P 2 — q 2 
plane. Finally, with regard to electron scattering, in general the singularity structure of the 
QPV in the time-like region of photon momentum squared dominates the behavior in the 
low momentum region of space-like q 2 . This feature is exhibited here despite the absence of 
a bound-state pole. The intercept at q 2 = 0, is provided by the constraints, (|16|) , derived 
from the WI, which the coefficient functions given in (|27D satisfy explicitly. Further out 
in the space-like region a transition to asymptotic behavior is expected. As is illustrated 
in Fig. 3b, an unphysically sharp transition occurs here due to the nonanalytic structure of 
the self-energy amplitude, B, in (^). Based on these observations, the approximation fl24|) 
does not provide a useful description of bound states or the approach to asymptopia, but is 
of pedagogical value for investigating the behavior of the QPV for the situation where the 
confinement length scale, defined here as the quark separation at which the vacuum becomes 
apparent, is taken to infinity. 



IV. NUMERICAL EVALUATION 

A. Methods and approximations 

The set of equations for the five transverse coefficient functions represented in (|17D , and 
given explicitly in the appendix, forms the basis for the numerical investigation undertaken 
here. The approach used in evaluating these equations is to express each element of the 
matrices (vectors) in the five-component space of coefficient functions as a matrix(vector) in 
the direct-product space of the magnitude of the momentum P and the direction cosine 
Cp q . The value of the photon momentum, q 2 , is a parameter of the equations. That is, 
for each value of q 2 a different set of equations is obtained. For the numerical evaluation, 
the momentum P and the direction cosine Cp q are discretized in terms of Gauss-quadrature 
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points. The integrations are then carried out as matrix multiplication. The coupled equa- 
tions are solved as a matrix equation in the expanded direct-product space of the coefficient 
functions and the Gauss integration points. The dimension of the equation is given by 
5xnpptsxncpts, where nppts and ncpts are the number of Gauss points representing the P 
and Cp q integrations respectively. The number of points in both integrations is varied to 
ensure accuracy. For the momentum, twenty to forty quadrature points are used, while six 
to ten are used for the direction cosine. 



1. Gluon two-point function ansatz 

The aim in this initial investigation is not necessarily to reproduce the known bound-state 
spectrum, but rather to correlate trends in the spectrum with changes in the parameteriza- 
tion of the infrared form of the gluon two-point function. The sensitivity of the spectrum 
to these changes then provides a constraint on the calculation of the QPV in the space-like 
region. To this end, a sufficient parametrization of the gluon two-point function is given by 

9„/„9n 37r 4 77 2 . R 2p2 M 167T 2 . . 

9 DiP » = ^T^- R ' P /4 + 11P n n{e + P y Ar < 28 > 

where the first term (which has the coordinate space representation e~ r2 / R °) provides a 
length scale, Rq, characterizing the infrared behavior while the second term models the 
known ultraviolet form. The parameter A, which is associated with the QCD scale, is 
fixed at 1 fm _1 throughout this investigation. As Rq approaches infinity the first term 
reduces to the delta-function model of Eq. (|24]) . The gluon two-point function fl28|) has been 
employed previously in the study of Bethe-Salpeter and Schwinger- Dyson equations || , and 
is expected to provide an acceptable accounting of the infrared effects for the purposes of 
the investigation conducted here. 
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2. Quark self energy ansatz 



The solution of the Bethe-Salpeter equation in the time-like region requires knowledge 
of the quark Green's functions, and hence the self-energy functions, for complex values of 
their arguments. The solution of the Schwinger- Dyson equation in the complex plane is a 
challenging problem in itself [13,7], and is being pursued elsewhere for forms such as that 
in (p8| ) Jljfl . Although these solutions are crucial to a comprehensive study of the present 
problem, reasonable approximations are available. The analytic solutions given in ( |25|) 
for the self-energy functions obtained using the simple model two-point function (^) are 
useful in this regard. It is shown, for example, in Ref. that the solutions to (|8|) and (|9|) 
for the self-energy amplitudes obtained using (|28|), display behavior similar to the simple 
model solutions of ( p5|) on the real axis with the exception of an additional tail on the scalar 
amplitude, B, attributed to the ultraviolet contribution in (^8|) . Others have argued on both 
theoretical and experimental grounds that the asymptotic form of the scalar contribution 
to the self energy should behave as B(P 2 ) ~ Am^/P 2 , which, in the Landau gauge, is exact 



within a logarithm |15|| . The self-energy functions used in the numerical work here are thus 
modeled as 

( 2, 
l + (l + 



A{P 2 



p2 



B(P 2 



p2 ) 



P 2 < \ 
— o 

p 2 > £ 



(29) 



where the quantity nip, has been referred to previously as the "dynamical-quark mass" 
P~5| , p76[1 , which sets the scale for dynamical chiral symmetry breaking. Its value is determined 
here by demanding that the function B and its first derivative match at the point rj 2 /6, and 



is thus given by mo = [24\/3]" 



1/3 



V 



0.29i]. For i] ~ 1 GeV, a dynamical-quark mass 



of mo ~ 290 MeV is obtained, which is consistent with the notion of a constituent-quark 
mass. In Ref. [[TJJ, using evidence from e + e~ annihilation, the dynamical mass is estimated 
to be mo k 244 MeV. These arguments place constraints on the range of the parameter 7], 
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which is thereby fixed at 5 fm _1 (« 1 GeV) for the remainder of this work. The ansatz ( |29|) 
encompasses the expected behavior of the solutions to the Schwinger-Dyson equation (|2|) on 
the real axis |7| , and is here taken to define the solutions in the portion of the complex plane 
sampled by the QPV calculation. 



B. Results 

1. Time-like q 2 and bound states 

Shown in Fig.[| are the solutions to Eq. (fT^) for the five transverse coefficient functions in 
the time-like region of photon momentum squared. Two different values of the infrared scale, 
Rq, are reported, from which it is observed that the spectrum is quite sensitive to changes 
in this parameter. In particular, an increase in i?o leads to an increase in the density of 
bound-state poles. This situation is, for example, analogous to a potential well, where Rq 
plays the role of the range while the barrier is provided by the repulsive interaction with the 
vacuum. With the value of R set to 1 fm, the ground-state mass occurs at 880 MeV. This is 
in reasonable agreement with the observed p mass, considering that pion dressing is expected 
to lower this value |DJ, yet the excitation spectrum does not coincide with experimentally 
observed quantities. Comparison with the measured excitation spectrum can be improved 
by decreasing the range, Rq, however this also leads to an increase in the ground-state mass, 



and away from the known value. As illustrated previously with the analytic solutions fl27|) , in 
the limit as Rq approaches infinity a continuous spectrum is obtained. From this discussion, 
and the behavior displayed in Fig.^j, it is evident that the combination of the gluon two- 



point function ( P8|) and the self-energy functions fl29|) , is not sufficient to reproduce the 
empirical excitation spectrum. Nevertheless, it is encouraging that this approach affords 
the investigation of spectra beyond the ground state, which until recently || has been the 
focus of Bethe-Salpeter phenomenology applied to mesons. The reproduction of spectra is 
clearly essential to constraining the interaction between quarks. In this regard, the utility 
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of the approach presented here is that, given a gluon two-point function and the associated 
solution to the Schwinger- Dyson equation, the vector bound-state spectrum can be obtained. 
The ability to reproduce the experimentally observed spectrum is then contingent upon the 
capacity of the ladder approximation to capture the relevant dynamics of QCD. 

The extraction of bound-state wavefunctions is achieved by considering the dependence 
of the coefficient functions, Xj , on the momentum P 2 and the direction cosine Cp q in the 
vicinity of a pole. Shown in Fig.^ are the solutions for the ground-state with mass 880 MeV, 
and with Rq = 1 fm, as a function of P 2 for the two limiting values of the direction cosine, 
Cp q = 0,1. The normalization is chosen such that the largest coefficient function is unity 
at P 2 = Cp q = 0. The relative strength of the other coefficients can then be observed. The 
absolute normalization is obtained through application of Eq.(^), but is not relevant to the 
present discussion. A rough estimate for the dependence of the solutions on the momentum 
and the direction cosine can be obtained by fitting the curves to an exponential of the form 
e~ p ( 1+aC W/ A ) which gives a ~ 0.29 and A ~ 3.5 fm^ 1 . The corresponding estimate of the 
spatial extent is ry ~ 2/A ~ 0.57 fm. 

2. Space-like q 2 and electron scattering 

The description of the QPV in the space-like region is largely independent of the spec- 
trum in the time-like region with the exception of the ground state, which dominates the 
low momentum structure. The approach to asymptotic behavior at higher momentum is 
primarily determined by the asymptotic forms of the gluon two-point function and the self- 
energy functions. Care has been taken here to accommodate both of these low and high 
momentum attributes in a manner that is consistent with their known features. The results 
of the calculation in the space-like region are summarized in Figs.[7] and |[ 

In Fig.0 the solutions for the coefficient functions, Xj , are displayed versus the photon 
momentum, q 2 , for P 2 = Cp q = 0. A smooth transition from the low momentum region, 
dominated by the ground state vector-meson pole, to the high momentum approach to 
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asymptopia is observed. It is therefore expected that the range of applicability of these 
solutions spans the entire range of momentum transfers available in electron scattering ex- 
periments. With the exception of A3 , which approaches unity asymptotically, the coefficient 
functions vanish for large momenta consistent with asymptotic freedom ||. The asymptotic 
form is apparent at q 2 ~ 50 fm~ 2 (~ 2 GeV 2 ), where the coefficient function A^ achieves 
roughly twice the strength of the other coefficients. It is also worth pointing out that to 
a very good approximation the coefficients A^ and A5 are equal over the momentum range 
reported. 

Shown in Fig.[8] are the coefficient functions versus the momentum P 2 for q 2 = 10 fm -2 , 
and for the limiting values of the direction cosine, Cp q = 0,1. Only a modest angular 
dependence is apparent, and again the coefficients Af and A5 are approximately equal. The 
observed behavior in the space-like region can be roughly characterized in terms of the scales 
of the problem by the following analytic forms: 

A - K " i+^i+^s.) I ln ( e + (p2 + « 2/4)/A2 )l " 2 l * 3 

A " K 1 + TTTJmI (e + (P '' + ?2/4,/A2) ) ~" 2 ' (30) 

where the constants Oj can be read from Fig.[7|, and in present case the ground state vector 
mass is M p = 880 MeV. The falloff with the momentum P 2 is reminiscent of the self-energy 
functions, while that with q 2 reflects the presence of the singularity in the time-like region. 
It should be understood that these forms are not the result of a detailed fit, but are rather 
offered as a guide to the qualitative behavior over the momentum range reported. 



V. SUMMARY 

The inhomogeneous Bethe-Salpeter equation in the ladder approximation has been used 
to study several issues that arise in the nonperturbative description of the electromagnetic 
interaction with quarks. The investigation presented here represents a first solution of 
this equation in a gauge invariant formulation which features confinement and asymptotic 
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freedom, in addition to the dynamical generation of qq vector meson modes. The principle 
results are summarized in Figs.|5|-^|. 

The solution in the time-like region warrants further study, particularly with regard to 
the reproduction of the excitation spectrum. Nevertheless, the description presented here 
shows promise in the ability to study the vector meson spectrum beyond the ground state. 
Further, the demonstrated sensitivity of the spectrum to the introduced infrared scale has 
positive implications for constraining the largely unknown behavior of the gluon two-point 
function at low momentum. The utility of this approach is, however contingent upon the 
ability of the ladder approximation to capture the relevant bound-state dynamics of QCD. 

The solution in the space-like region is potentially the most useful result of this investiga- 
tion. The quark-based study of the EM properties of hadrons relies heavily on the capacity 
to describe the quark-photon interaction. Modeling this quantity based on gauge invariance 
provides only partial constraints on its behavior, and in particular allows for arbitrariness 
in the components transverse to the four momentum of the photon, which are those probed 
by electron scattering experiments. The use of vector meson dominance or the bare QPV, is 
limited by the momentum range for which these approximations are sensible. Here a some- 
what more fundamental approach has been taken by describing the quark-photon interaction 
through a model in which the dynamics dictate the role of qq vector bound states and the 
asymptotic behavior. The transverse contributions to the QPV are thus uniquely defined 
within the employed ladder approximation to the inhomogeneous Bethe-Salpeter equation, 
and no restrictions on the momentum range are implied. The solutions obtained here are 
therefore expected to provide an essential ingredient in the quark-based description of the 
EM interactions with hadrons. 
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DETAILS OF EQ.@ 

The essential steps required to reduce the inhomogeneous Bethe-Salpeter equation flU) 
to the matrix form given in (|T7|) are listed below. To begin, the coefficient functions are 
isolated using ([12]) and (0), and the orthogonality of the Dirac matrices to obtain 



(pT\2 i pT 

^Ar = i^tr[rv(p, 5 )] 

i^Aj - 3X1 = ~ (V - "-f) tr [7J„(P, ,)] 

(pT\2 jP T P T 

~2 -^2 - AS 1 = -^^tr [ 7l/ r,(P, g)] (31) 
+ = -l^tr [ 7 ^(P,,)] 

—-4 — A 5 = --e^ap— 2-tr bWsiyP, g)] . 

7] O 

The right-hand side of the equations in ( plf ) are evaluated using the inhomogeneous Bethe- 
Salpeter equation (|1|) for r M . 

The remaining obstacle is the evaluation of the four- dimensional integration present in 
(H), which, upon substitution into fl3ll) , produces expressions of the form 

1=J d 4 KD (P 2 + K 2 -2PK^jl - C 2 Pq C KT -2PKC Pq C Kq ) 

xC n KT F(p 2 ,K 2 iq 2 ,C Kq ,C Pq ), (32) 

where n = 0,1,2; CKqiCxr) is the direction cosine between K and q(K and P T ); and 
D ((P — K) 2 ) = g 2 D ((P — K) 2 ) I (3r] 2 7T A ). Here F is a generic function representing the 
integrands obtained from (|3T|). Since the vectors Pj and g M are orthogonal, they can be 
used to define two of the four axes of integration. Rewriting the integration as 

//*oo rl 
d 4 K = 2 j dKK 3 J ^ dC Kq dC KT dC 1 dC 2 S{l - C 2 Kq - C 2 KT - C\ - C 2 ), (33) 

where C\ and C 2 are the remaining direction cosines, then obtains from ( |32"D the result 
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l = 2n J™ dKK 3 dC Kq F (P 2 , K\ q 2 , C Kq , C Pq ) D n (P 2 , K 2 , C Pq , C Kg ), (34) 

where 

D n (P 2 , K 2 , C Pq , C Kq ) = J* dC KT D (P 2 +K 2 -2PK^l-C 2 Pq C KT -2PKC Pq C Kq ) C n KT) 



(35) 



and x = a/1 — C\ q . 



The explicit form of the quantities in (|T7j ) is listed below. The vector b is given by 

b= (0,3, 1,0,0). (36) 



The matrix M is given by 



Mi = l A4 2 = -^ 



AT 23 = 3 AT 32 = -^ 
A/" 33 = 1 A/42 = 1 

A/44 = A/55 — 1) 



where all other elements are zero. The matrix M. is given by 

M n = InDiFofr M 12 = 2-kD^ (v + $C 2 Kq W) 

M 13 = -2-nb lp \V M u = 2nD 1 ^LC 2 Kq (V + 



M 21 = -7rD fx 2 V M 22 = nbo^x 2 (Fx - 2^T 

M 23 = -nb (3F X - 2^T) M 24 = -2nD ^fix 2 C 2 Kq T 

M 25 = 2nt) ^fx 2 T M 31 = -nD 2 ^V 
M 32 = nD 2 ^ (F x - 2^T) = -7T (F X D Q - 2^TD 2 

A4 34 = -2nD 2 ^fc 2 Kq T M 35 = vr (b x 2 - b 2 ) ^fr 

M 41 = 7ib 17 ^C Kq (v + $w) M i2 = -^b x ^rC Kq F 2 

A4 43 = -2TX b l p Lp T C Kq T M44 = —TvDi pLp T CKq^F 3 

M 53 = -nD^T M 55 = nb^Fo 
where all unlisted elements are zero, and 



(37) 



(38) 
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K 2 - q 4 ) 4t + u 



K 2 Cl q 



31 
■l 



(39) 



F = - K') ±T + U F 
F 2 = F 1 + 2 (£ - F 2 ) 4, F 3 = Fx + 2 (s 

T = r/ 4 a(F 2 )a(Fi) C/ = 7] 2 (3(K 2 + )(3(K 2 ) 

V = rf [a(K+)P(K-) + a(F_)/3(F+)] W = rf 

The quantities a and (3 are defined in ([XO|) ; T,U,V, and W are dimensionless functions of 
K 2 ,q 2 and Cj^ ; and the components of the momentum F are defined as Fj, = PCp q and 
F T 



g(^ + )/3(^-)- a (A--)/3(^ + ) 
2Kq 



PJI 



r<2 

Lyp q . 
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FIGURES 

FIG. 1. The EM vertex for a composite pion at tree level. The quark Green's functions, QPV, 
and the pion Bethe-Salpeter amplitudes, are dressed in a consistent, gauge-invariance manner. 

FIG. 2. The ladder approximation for (a) the QPV and (b) the quark self energy. The Green's 
functions are defined by their inverse, G~ 1 (P) = i P + S(P). 

FIG. 3. The solutions for the five transverse coefficient functions, Xj , obtained using the simple 
gluon two-point function described by Eq,(|2~4|) in the text, are plotted as a function of the photon 
momentum q 2 in both the (a) time-like and (b) space-like regions for P 2 = Cp q = 0. 

FIG. 4. The coefficient function \^ is plotted in the q 2 — P 2 plane for Cp q = 0. The variation 
of the singularity with P 2 is indicative of a continuous spectrum. 

FIG. 5. The five transverse coefficient functions are plotted versus time- like photon momentum 
q 2 for P 2 = Cp q = 0. The dependence upon the infrared scale, Rq, is illustrated by comparing (a) 
where Rq = 1 fm with (b) where Ro = 0.75 fm. 

FIG. 6. The five transverse coefficient functions are plotted versus the momentum P 2 for values 
of the direction cosine (a) C Pq = and (b) C Pq = 1, and for -q 2 = My = (880 MeV) 2 . 

FIG. 7. The five transverse coefficient functions are plotted versus the photon momentum q 2 
for P 2 = C Pq = 

FIG. 8. The five transverse coefficient functions are plotted versus the momentum P 2 for photon 
momentum q 2 = 10 fm -2 , and the limiting values of the direction cosine (a) Cp q = and (b) 
Cp q = 1. 
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